Phase transitions induced by variation of invasion rates 
in spatial cyclic predator-prey models with four or six species 
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Cyclic predator-prey models with four or six species are studied on a square lattice when the invasion rates 
are varied. It is found that the cyclic invasions maintain a self-organizing pattern as long as the deviation of 
the invasion rate(s) from a uniform value does not exceed a threshold value. For larger deviations the system 
exhibits a continuous phase transition into a frozen distribution of odd (or even) label species. 



PACS numbers: 87.23.Cc, 89.75.Fb, 05.50.-l-q 
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Spatial predator-prey models with cyclic competition ex- 
hibit several peculiar features that can be relevant in biologi- 
cal, chemical and ecological systems d Q, B 0. Btl . In these 
systems species are distributed on the sites of a lattice and 
the predators' offspring can be substituted for the neighbor- 
ing prey. More precisely, for the simplest cyclic n-species 
systems (n > 3) a randomly chosen individual of species 
i can occupy the neighboring site of species i + 1 mod n 
(i = 0,l,...,n — 1) with an invasion rate independent of i. 
On a one-dimensional lattice (d = 1) this evolutionary pro- 
cess yields growing domains (of sites occupied by the same 
species) if n < 5; otherwise the system develops into a frozen 
pattern |6]. For higher-dimensional lattices the uniform inva- 
sion rate maintains a self-organizing pattern if the system size 
N goes to infinity and the number of species does not exceed 
a threshold value dependent on d, that is, n < nc{d) [e.g., 
nc{d = 2) = 14] [7J. In this case the self-organization is 
maintained by the unceasing motion of interfaces separating 
a predator from its prey and the species are present with the 
same average density {pi = 1/ n) independently of the initial 
state containing all species. 

The mean-field analysis of these systems indicates differ- 
ent types of solutions [|4i]. Besides the above mentioned sym- 
metric solution (pi = 1/n) these models have n homoge- 
neous solutions (e.g., po = 1 and pi = ... = p„_i — 0) 
which are unstable against the invasion of the corresponding 
predator. At the same time this system exhibits a continuous 
set of oscillatory solutions satisfying conservation laws, e.g., 
J2iPi — 1 ^^'^ YiiPi — constant |4]. The pair approxima- 
tion, as a first step to extend the classical mean-field theory 
IllH], confirms the existence of both the symmetric and the n 
homogeneous solutions. Instead of the periodic oscillations, 
however, the pair approximation predicts growing oscillations 
ending in one of the homogeneous states due to the numerical 
rounding or noise. For the simplest case, called the evolu- 
tionary rock-scissors-paper game (n — 3), the latter behavior 
was observed on Bethe lattices of degree 2; > 6 in the limit 
N ^ 00 ^ [l^l (while a limit cycle is reported for z — 3 
or 4). The real behavior, observed by means of Monte Carlo 
(MC) simulations, on a square lattice (for n — 3) could be 
reproduced by the more sophisticated four-site approximation 
determining the probability of each possible configuration on 
the 2x2 cluster of sites ^ The failure of the pair ap- 
proximation for spatial models indicates that their behavior is 
sensitive to the details of the lattice structure. Similar relevant 



sensitivity is found when considering the effect of different 
(species-specific) invasion rates. 

When the simplest (n = 3) case is considered, the varia- 
tion of invasion rates modifies the composition of the central 
(stationary) solution |4]. Paradoxically, if one of the species 
is fed more effectively by increase of its invasion rate then its 
predator will benefit from this modification | lll [T2ll . In other 
words, if the largest invasion rate is from species to 1 then 
species 2 (the predator of 0) will have the highest density in 
the central stationary state. This mean-field prediction is con- 
firmed by MC simulations on two-dimensional lattices. The 
unusual response of the system is related to the indirect in- 
teraction mediated through the food web. For example, if the 
density of species is increased by any external effect then it 
yields a decrease in the density of species 1, which consume 
less of species 2, and so on. For a cyclic food web with odd 
number of species this mechanism results in a negative feed- 
back, in agreement with the mentioned results in the station- 
ary state. For even number of species, however, the feedback 
becomes positive and may yield the extinction of the odd (or 
even) label species. 

Using mean-field analysis, pair approximation, and MC 
simulations, this parity effect was demonstrated by Sato et al. 
ifTsIl who considered a cyclic predator-prey model where each 
invasion rate is chosen to be unity excepting the rate {a) from 
species to 1 . According to their results for even n the odd 
label species die out if a < 1, while the even label species 
become extinct if a > 1. Here a = 1 denotes a singular point 
where the symmetric solution (pi ~ l/n) remains stable. Un- 
fortunately, the authors did not study the close vicinity of this 
singular point in their MC simulations. 

Now we show that, instead of sudden changes in pi at 
a = I, the systems exhibit two consecutive phase transitions 
referring to the existence of an n-species state with continu- 
ously varying composition. We also study the critical behav- 
ior of the phase transition separating the frozen and coexisting 
states. 

In order to improve accuracy, first we consider the sym- 
metric version of the model suggested by Sato et al. for 
n — 4. It should be mentioned that the qualitative results 
obtained for the symmetric model agree with those valid for 
the original model as discussed later. We assume that each 
site a; of a square lattice is occupied by a single individual 
belonging to one of the n species and their spatial distribu- 
tion is described by a set of site variables Sx = i where 
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i = 0,...,(n — 1) refers to the label of the resident species. 
The system evolution is governed by the repetition of invasion 
between two neighboring sites chosen at random if the sites 
are occupied by a predator-prey pair [i.e., {s^, Sy) — > {s^, s^) 
if Sx is the predator of Sj,] . The species-specific invasion prob- 
abilities are proportional to the invasion rates denoted in Fig.[T| 
and backward invasions are forbidden. Notice that if a suit- 
able new time scale is chosen for a > 1 this system becomes 
equivalent to a model with a < 1 if the species labels are 
increased by 1. 




FIG. 1 : Invasion rates between the predator-prey pairs for the four- 
species cyclic food web. 

MC simulations are performed for different values of a on 
a square lattice with periodic boundaryy conditions and the 
linear size L of the system is varied from 600 to 4000. The 
large system size is chosen to avoid artifact effects like trap- 
ping to an absorbing state as a consequence of small system 
size. During one time unit, called a MC step (MCS), each site 
has a chance once on average to be invaded and to invade one 
of the neighboring sites. The system is started from a random 
initial state and after a suitable thermalization time tt the av- 
erage density of species (pi) is determined by averaging over 
a sampling time (ts). The thermalization and sampling times 
(tt ~ ^ 5Tr where is the relaxation time) are varied from 
10 000 to 200 000 MCSs. Clearly, the largest values of L, U, 
and ts are used in the close vicinity of the transition points 
where the relaxation time, correlation length, and fluctuations 
diverge algebraically. Due to the translational symmetry in 
the invasion rates, the odd (even) label species have the same 
density, that is, po = p2 and pi = ps = 1/2 — pQ. Henceforth 
our analysis will be focused on the limit N ^ oo. 

According to the MC results (see Fig. |2|, two consecutive 
critical transitions can be observed for = 0.9332(1) and 
= 1/ «ci ■ For a < ad the simulations end in a frozen 
spatial distribution of species 1 and 3. Clearly, the even la- 
bel species (0 and 2) form a frozen pattern in the final state if 
a > ■ Both types of frozen pattern can be considered as 
an infinitely degenerate absorbing state from which the system 
cannot escape. Within the intermediate region the four species 
coexist and invade each other cyclically. In this phase, the av- 
erage density of species varies monotonically with a as illus- 
trated in Fig. 121 This feature is well reproduced by generalized 
mean-field techniques where we have determined the proba- 
bility of each possible configuration on 2 x 2- or 3 x 3-site clus- 
ters on a square lattice (for a brief description of this method 
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FIG. 2: Average density of species as a function of a characterizing 
the invasion rates from the even label species as denoted in Fig. [T] 
Diamonds and squares indicate MC results while the dashed line 
shows the prediction of the pair approximation [13] for the even la- 
bel species. The results of 2 x 2- and 3 x 3-site approximations are 
illustrated by dash-dotted and solid lines. 

see Appendix C in Ref. fs*]). Clearly, the analytical approxi- 
mations predict mean-field type transitions and the numerical 
values for the first critical point are: a^t^' — 0.9267(2) and 
ag^^ = 0.9252(1). 

The rigorous analysis of the MC data in the vicinity of 
the first transition point indicates critical behavior Namely, 
the density of the vanishing species, as an order parameter, 
follows a power law behavior when approaching the critical 
point, that is, po — p2 — {a ~ Uc^Y where (3 = 0.58(1) as 
illustrated in Fig. |3] Similar behavior characterizes systems 
belonging to the directed percolation (DP) universality class. 
In fact, when the transition point is approached, the extinction 
of the four-species cyclic phase for sufficiently large time and 
length scales becomes analogous to the elimination of infected 
sites in the contact process 1 14]. 

The universal behavior of the directed percolation uni- 
versality class ifTsIl occurs typically on homogeneous back- 
grounds |8]. In the present case, however, the extinction pro- 
cess leaves an inhomogeneous, frozen pattern behind. Conse- 
quently this system has infinitely many absorbing states. Con- 
sidering a much simpler system, Jensen 1 16] has demonstrated 
that in such systems the non-equilibrium phase transition may 
also belong to the directed percolation universahty class. 

The present MC results give another example where the 
universal behavior is not affected by some types of inhomo- 
geneities. In other words, the extinction process is controlled 
by a background that can be considered homogeneous on suf- 
ficiently large time and length scales. To confirm this state- 
ment, we have also determined the time dependence of the 
density of the vanishing species at the first critical point. It is 
well known that for the two-dimensional, directed percolation 
universality class the order parameter vanishes algebraically, 
that is, pQ{t) = p2{t) ~ where 9 = 0.45(1) |17]. The 
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FIG. 3: Log-log plot of the vanishing density of species vs. a — aci ■ 
The dashed line with a slope of 0.58 indicates the power law behavior 
characteristic of the directed percolation universality class in two- 
dimensional systems. 



MC data (see Fig. justify this expectation. 
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FIG. 4: Log-log plot for the time dependence of the density of the 
vanishing species at the critical point (L = 3000). Fluctuations in 
po (t) — p2 (t) are suppressed by averaging over a time window with 
a typical width of t/6. The dashed line (with a slope of -0.45) illus- 
trates the power law behavior characteristic of the directed percola- 
tion universality class. 

The existence of the four-species phase (and phase transi- 
tions) is independent of the symmetrical invasion rates illus- 
trated in Fig. [T] According to the MC simulations the four- 
species phase is also observed for those models where only 



one invasion rate differs from 1. Following the notation of 
Sato et al. [i3J, we have studied the system where the inva- 
sion rate from species to 1 (denoted by a) is varied while the 
other rates remain unity. Here, in the absence of the previously 
mentioned symmetry all the species densities become differ- 
ent within the four-species phase existing if < a < ctca 
[here = 0.8716(1) and Uc^ = 1.145(1)]. When is 
approached from above, the densities of even label species 
decrease simultaneously, i.e., po, p2 — {a — a^)^ (po > P2), 
suggesting that the critical behavior belongs to the DP uni- 
versality class too. An asymmetric composition can also be 
observed in two-strategy frozen patterns, namely, pi > p3 if 
a < , and p2 > po if a > . Similar behavior is ex- 
pected when all the invasion rates are tuned. The robustness of 
this universal behavior was confirmed for those systems where 
a slight mixing is introduced by allowing site exchange with 
a low probability between neighboring neutral species (e.g., 
and 2 or 1 and 3) Sill. 

In addition to the four-species systems, we have also stud- 
ied what happens in a six-species cyclic predator-prey model 
if the invasion rate a from species to 1 is varied while the 
others are chosen to be unity. In the spirit of the parity law 
ifisll . the cyclic food web mediates an indirect effect yield- 
ing the behavior described above. MC simulations have con- 
firmed that this system evolves from the random initial state 
to a frozen pattern consisting of odd (even) label species if 
a < aci = 0.8971(3) [a > a^^ = 1.1196(3)]. Within the 
intermediate phase (aci < a < ) all six species survive 
and their densities vary similarly to those plotted in Fig. |2] 
According to our previous paper [20], similar behavior is ex- 
pected for n = 6 if some weak mixing is allowed, while the 
unusual behavior within a range of the strength of mixing for 
n — and 10 raises further questions. 

In summary, the present investigation revises the parity law 
deduced originally from the results of the pair approxima- 
tion for spatial cyclic predator-prey models with even number 
of species. The more sophisticated levels of the dynamical 
mean-field approximation and also the MC simulations show 
that a (noisy) cyclic invasion is capable of maintaining a state 
in which all the species survive with a composition depend- 
ing on the invasion rates if the unilateral deviation from the 
uniform rate does not exceed a threshold value. When the 
invasion rate(s) is varied the system exhibits (directed perco- 
lation type) critical phase transitions into frozen states where 
only the odd (or even) label species are present. 
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